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Abstract 

A  large  set  ofquasi-classical,  direct  dynamics  trajectory  simulations  w  ere  performed  for  decomposition  of 
1,  5-dinitrobiurct  (DNB)  over  a  temperature  range  from  4000  to  6000  K.  aimed  at  providing  insight  into 
DNB  decomposition  mechanisms.  ITtc  trajectories  revealed  various  decomposition  paths,  and  reproduced 
the  products  (including  HNCO,  NsO,  NO>,  NO.  and  water)  observed  in  DNB  pyrolysis  experiments. 

Using  trajcctoiy  results  as  a  guide,  stmcturcs  of  intermediate  complexes  and  transition  states  that  might  be 
important  for  decomposition  were  determined  using  density  function  theory  calculations.  Rice- 
Ramspcrgcr-Kasscl-Marcus  (RRKM)  theory  was  then  utilized  to  examine  behaviors  of  the  energized 
reactant  and  intermediates  and  to  determine  unimolccular  rates  for  crossing  the  various  transition  states. 
According  to  RRKM  predictions,  the  dominant  initial  decomposition  path  of  energized  DNB  corresponds 
to  elimination  of  HNNO:H  via  a  concerted  mechanism  where  the  molecular  decomposition  is 
accompanied  with  intra-molecular  H-atom  transfer  from  the  central  nitrogen  to  the  terminal  nitro  oxygen. 
Other  important  paths  correspond  to  elimination  of  NO:  and  H;NNO:.  NO:  elimination  is  a  simple  N-N 
bond  scission  process.  Formation  and  elimination  of  nitramidc  is  however  dynamically  complicated. 
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requiring  twisting  a  -NHNO2  group  out  of  the  molecular  plane,  followed  by  an  intramolecular  reaction  to 
form  mtramidc  before  its  elimination.  These  two  paths  become  significant  at  temperature  above  1500  K. 
accounting  for  >  17%  of  DNB  decomposition  at  2000  K.  This  work  demonstrates  that  quasi-classical 
trajectory  simulations,  in  conjunction  with  electronic  structure  and  RRKM  calculations,  are  able  to  extract 
mechanisms,  kinetics,  dynamics  and  product  branching  ratios  for  the  decomposition  of  complex  energetic 
molecules,  and  predict  how  they  vary  with  decomposition  temperature. 
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I.  Introduction 


1,  5-dinitrobiurct  (DNB,  C2H3N5O6,  molecular  weight  193 .OX)  is  a  very  powerful  and  promising  new 
explosive.1*4  This  nitrogen-rich  compound  was  first  synthesized  by  Thiele  and  Uhlfcldcr  in  I898.5 
Recently.  Klapotkc  and  co-workers  reported  both  experimental  and  theoretical  studies  of  DNB,  which 
determined  DNB  molecular  structure,  energetics,  pyrolysis  as  well  as  detonation  velocity  and  pressure.1 
DNB  has  an  approximate  planar  structure,  and  the  standard  heat  of  formation  (ArH)  was  calculated  to  be 
-353  IJ/mol  (-1,829  kJ/kg)  at  the  MP2/cc-pVTZ  level  of  theory.2  The  enthalpy  (A,H)  for  the  reaction  of 
DNB(s)  ->  2CO:(g)  +  3/2  H:(>(l)  +  1/4  02  (g>  +  5/2  Nj(g)  was  determined  to  be  1,003  ±  39  kJ’mol 
(5.195  ±  200  kJ/kg)  using  bomb  calorimetry.2  DNB  detonated  strongly  in  the  steel  sleeve  test,  without  the 
need  of  adding  any  additional  oxidizer.  The  detonation  velocity  and  detonation  pressure  were  calculated 
to  be  8.66  mm(is  and  33.9  GPa,  respectively.*'*  These  values  arc  close  to  those  of  the  well-established 
energetic  materials  such  as  HMXh  (i.c.  cyclotctramcthylcne-tctranitraminc)  and  RDX6  (i.e. 
cyclotnmethylenc-trinitramine).  Thermal  decomposition  of  DNB  was  examined  by  themiogravimctric 
analysis  (TGA)  and  differential  scanning  calorimetry  (DSC),  and  gaseous  products  were  identified  using 
mass  spectrometry  and  Fourier  transform  infrared  spectroscopy.1 4  DSC  analysis  revealed  that  thermal 
decomposition  of  DNB  docs  not  occur  in  a  single  step:  instead  it  has  several  ditTcrcnt  activated  reactions. 
The  ultimate  decomposition  gases  include  UNCO.  N:(),  CO2,  and  water.  However,  the  decomposition 
mechanism  and  dynamics  have  not  yet  been  clearly  established.  One  of  the  challenges  is  that  the 
mechanisms  involving  highly  reactive  intermediates  arc  difficult  to  probe  experimentally,  since  such 
dissociation  is  a  violent  event  and  over  within  nanoseconds. 

It  is  worth  noting  that  Chambreau  et  a/,  reported  a  hypcrgolic  ionic  liquid  fuel  based  on  dicyanamidc 
anions  with  white  fuming  nitric  acid.  Hypcrgolic  systems  arc  fuel-oxidizer  combinations  that  ignite 
spontaneously  upon  mixing  of  two  propellants  at  low  temperatures  and  pressures,  without  the  presence  of 
an  external  ignition  source.  According  to  the  proposed  mechanism  for  the  hypcrgolic  ignition  of 
dicyanamidc  anions  with  nitric  acid,  the  reaction  sequence  results  in  the  formation  of  DNB  as  an 
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intermediate  in  the  ignition  process.**  Therefore,  a  detailed  understanding  of  mechanisms  and  dynamics  of 
DNB  decomposition  is  necessary  for  constructing  an  accurate  kinetics  model  to  predict  the  fuel 
performance,  for  optimizing  the  design  of  propulsion  systems  using  hypergolic  bipropellants,  and  for 
modifying  and  designing  new  fuels. 

In  the  present  study,  thermal  decomposition  of  DNB  was  probed  theoretically.  For  complex 
molecules,  especially  tor  those  with  high  internal  energy  content,  use  of  chemical  intuition  to  predict 
reaction  pathways  can  prove  unreliable.  There  may  be  concerted  reactions  that  are  difficult  to  predict, 
and  once  energy  stored  in  the  molecule  begins  to  release,  the  system  docs  not  necessarily  follow  the 
minimum  energy  reaction  path.  i.c..  the  subsequent  behavior  is  controlled  by  dynamics/  Even  for  steps 
in  the  reaction  that  arc  controlled  by  statistical  factors,  standard  transition  state  locating  methods  may  be 
difficult  to  apply,  because  the  molecular  transformations  are  too  complex.  A  useful  approach  to  treating 
such  system  is  quasi-classical.  direct  dynamics  trajectory  simulations,  where  the  motion  of  the  molecule  is 
followed  allowing  the  molecule  to  "show  us"  what  the  preferred  reaction  pathways  arc.  The  direct 
dynamics  method  dispenses  with  the  potential  energy  surface.  Instead,  it  calculates  the  energies,  force 
constants,  and  Hessian  "on  the  fly"  using  quantum  chemistry  methods/  This  method  becomes 
computationally  attractive  when  the  dimensionality  of  the  system  increases,  particularly  for  the  DNB 
molecule  which  contains  13  heavy  atoms.  Dynamics  simulations  partition  the  energy  generated  by 
exothermic  reactions  to  vibrational,  rotational  and  translational  modes,  increasing  the  chance  of  locating 
new  reaction  pathways.  In  addition,  by  following  the  variation  of  the  potential  energy  during  the 
trajectory  rather  than  relying  on  intuition,  we  can  identify  better  geometries  for  I  S  searching,  as 
demonstrated  below. 

The  remaining  of  the  manuscript  is  organized  as  follows.  Section  II  describes  the  methods  of 
electronic  structure-  statistical  calculations  and  quasi-classical  trajectory  simulations  used  to  determine  the 
thermal  decomposition  features  of  DNB.  Trajectory  simulation  results  arc  presented  in  section  III, 
including  an  overview  of  the  nature  and  time  scales  of  trajectories  we  arc  simulating,  the  prescription  for 
sorting  trajectories,  and  contributions  of  various  decomposition  mechanisms  at  different  temperatures. 
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Statistical  modeling  results  arc  presented  in  section  IV,  starting  with  construction  of  the  reaction 
coordinate  for  DNB  decomposition  using  density  function  theory  electronic  structure  calculations, 
proceeding  to  decomposition  product  branching  using  Ricc-Ramspcrgcr-Kasscl -Marcus  (RRKM) 
theory.1’  and  finally  to  the  analysis  of  inter-conversions  of  primary  decomposition  products.  Conclusions 
arc  presented  in  section  V. 

II.  Computational  methodology 

A.  Electronic  structure  and  RRKM  calculations 

To  aid  in  reaction  coordinate  interpretation  and  get  energetics  information,  density  function  theory 
(DFT)  electronic  structure  calculations  were  performed  at  the  B3 LY P/6-3 1++G**  level  of  theory,  using 
the  (iaussian  09  (B.OI )  suite  of  programs. Iv  Geometries  were  optimized  by  calculating  force  constants  at 
every  step.  Vibrational  frequencies  and  zero-point  energies  (ZPE)  were  scaled  by  a  factor  of 0.955  and 
0.9X1. ,J  respectively.  The  corrected  ZPE  were  added  to  the  stationary  point  energies.  All  the  transition 
states  (TSs)  found  were  verified  to  be  first-order  saddle  points  by  frequency  calculations,  and  the 
vibrational  mode  with  the  imaginary  frequency  to  correspond  to  the  reaction  coordinate.  When  necessary 
intrinsic  reaction  coordinate  (IRC)  calculations  were  used  to  determine  which  minima  are  connected  by  a 
particular  TS.  To  obtain  more  accurate  energies  of  the  reactant  and  products,  calculations  were  also 
performed  at  the  MP2/cc-pVDZ  level  of  theory. 

A  transit  ion- state- theory -based  statistical  model.  Ricc-Ranisperagcr-Kassel-Marcus  (RRKM) 

theory.1'  was  employed  to  examine  the  properties  of  the  reactant,  intermediate  complex  and  transition 

states  as  a  function  of  the  total  internal  energy  and  associated  angular  momentum.  The  fundamental 

assumption  in  RRKM  model  is  that  energy  is  randomized  and  distributed  statistically  among  all  the 

energetically  accessible  states  of  the  system,  and  the  rate  of  a  particular  process  (e.g..  dissociation  of  a 

molecule,  or  inter-conversion  between  different  isomers)  is  proportional  to  the  total  number  of 

energetically  accessible  states  at  the  transition  state.*1  As  a  consequence,  statistical  reactions  tend  to 

occur  by  paths  close  to  the  minimum  energy  path,  as  the  density  of  states  is  highest  for  such  paths.*' 

RRKM  unimolecular  rates  for  crossing  various  transition  states  and  density  of  states  were  calculated  with 
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the  program  of  Zhu  and  Hasc,24  using  its  direct  count  algorithm,  and  scaled  frequencies  and  energetics 
from  the  DFT  calculations.  The  angular  momentum  quantum  number./  used  for  RRKM  calculations  was 
determined  from  the  corresponding  DNB  rotational  energy,  assuming  DNB  can  be  roughly  treated  as  a 
nearly  symmetric  top  molecule  since  two  of  the  principle  moments  of  inertia  arc  not  very  different. 
Product  branching  is  given  by  the  ratio  of  the  RRKM  rates. 

B.  Quasi-classical  direct  dynamics  trajectory  simulations 

In  this  study,  direct  dynamics  trajectory  simulations  were  conducted  for  two  purposes:  1 ) 
investigating  decomposition  mechanisms  of  DNB  at  various  temperatures,  and  discovering  decomposition 
products;*'  and  2)  locating  key  transition  states  for  important  decomposition  paths,  when  the  standard  TS- 
scarching  methods  (i.c.,  TS.  QST2,  and  QST3  methods  in  Gaussian)  fail.**  *  Trajectory  simulations 
started  at  the  equilibrium  geometry  of  DNB.  A  general  chemical  dynamics  program  VENUS  of  Hasc  el 
ai  *K  was  used  to  set  up  the  trajectory  initial  conditions,  and  the  Hessian-based  method  of  Bakkcn  et  at. 
implemented  in  Gaussian  was  used  to  propagate  each  trajectory  with  Hessians  recalculated  every  five 
steps.  'Hie  integrators  were  performed  with  a  step  size  of  -  0.4  fscc.  which  conserv  ed  the  total  energy  to 
better  than  10*4  Hartrcc.  The  SCF  =  XQC  option  was  used  during  trajectory  integration  so  that  a 
quadratically  convergent  Hartrcc-Fock  (QC-SCF)  method1  '**  was  adopted  in  case  that  the  usual  but  much 
faster,  first-order  SCF  method  did  not  converge  within  the  allotted  number  of  cycles.  Because  millions  of 
gradient  and  Hessian  evaluations  were  required,  the  level  of  theory  used  for  simulations  was  necessarily 
modest.  To  select  a  suitable  level  of  theory,  we  calculated  the  enthalpy  of  combustion  of  DNB,*  i.c..  A,H 


for  C;H,N5(>6<s)  ->  2CO:(g)  +  3/2  H:0(l)  +  1.4  02  (g)  +  5/2  N:(g).  using  B3LYP  and  MP2  methods  with 
6-3 1G.  6-3 1G*.  6-3 1+G*.  and  6-31++G**  basis  sets.  The  enthalpy  of  vaporization  of  water  (44.00 
kJ/mol)  was  taken  from  the  NIS  I  data  base.’0  and  the  enthalpy  of  sublimation  of  DNB  was  assumed  to  be 
close  to  that  of  biuret2  (which  is  -126.7  kJ/mol  from  NIST30).  We  then  compared  these  results  to  the 
reported  benchmark  calculation  at  electron  correlated  MP2/cc-pVTZ  level  of  theory.*  and  the  results  are 
listed  in  fable  1.  It  was  found  that  the  MP2  method  ran  into  a  convergence  problem  for  DNB  with  the 
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selected  basis  sets.  Compared  to  MP2/cc-pVTZ.  all  B3LYP  methods  underestimate  A,H;  however, 
including  polarization  in  basis  sets  significantly  improves  the  accuracy  of  DFT  energetics.  (Note  that 
none  of  the  calculated  A,H  values,  including  the  MP2/cc-pVTZ  value,  perfectly  matches  the  experimental 
determined  enthalpy  of  detonation  of  1.003  kJ/mol.  which  is  most  likely  due  to  uncertainty  arising  from 
the  estimated  enthalpy  of  DNB  sublimation  and  different  structures  of  DNB  in  the  gas  phase  and  in  the 
crystalline  state.*)  On  the  basis  of  the  overall  level  of  agreement  and  computational  speed,  the  B3LYP/6- 
3  1G*  level  of  theory  was  chosen  for  the  mam  set  of  trajectories.  It  is  noted  that  the  major  purpose  of  our 
trajectory  simulations  is  to  probe  the  gross  features  of  DNB  decomposition,  and  to  identify  important 
decomposition  paths. 

The  initial  vibrational  and  rotational  energies  of  DNB  were  sampled  from  Boltzmann  distributions  at 
specified  temperatures.  The  quasi-classica!  initial  vibrational  states  were  simulated  by  giving  each 
reactant  atom  displacement  from  equilibrium  and  momentum  appropriate  to  the  initial  rovibrational  state, 
with  random  phases  for  the  different  modes.  The  molecule  has  zero-point  energy  in  all  vibrational 
models.  Because  the  decomposition  time  scale  at  typical  detonation  temperature  would  be  far  too  long 
for  practical  trajectory  simulations,  simulations  have  to  be  carried  out  at  high  temperature.  To  identify  a 
reasonable  temperature  range  for  dynamics  simulations,  we  first  ran  a  series  of  trajectories  at  a  low  level 
of  theory  HF/6-31G,  at  a  temperature  range  of  1000  K  to  8000  K.  Within  a  trajectory  simulation  period 
of  4  pscc.  no  dissociative  trajectories  were  observed  at  temperature  below  3000  K.  At  temperature  of 
3000  K.  only  15%  of  total  trajectories  were  dissociative.  On  the  other  hand,  almost  every  trajectory 
underwent  decomposition  within  a  short  time  (<  0.4  pscc)  at  temperatures  higher  than  6000  K.  Therefore, 
three  temperatures.  4000,  5000  and  6000  K.  were  chosen  for  running  trajectory  simulations  at  the 
B3LYP/6-31G*  level  of  theory.  For  each  temperature,  batches  of  trajectories  (50  trajectories  each)  w  ere 
calculated.  Trajectories  were  propagated  for  a  preset  length  of  time  (-  4  pscc).  or  until  the  distance 
between  the  decomposition  products  exceeded  9.0  A.  A  total  of  150  trajectories  were  completed  at  the 
B3LYP/6-31G*  level  of  theory  in  the  courec  of  the  study  reported  here,  each  taking  -  380  CPU  hours  on 
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an  Intel  Core  2  dual  and  Core  2  quad  (3.0  GHz)-bascd  64  bit  Linux  computational  cluster.  Thus  the 
complete  set  of  trajectories  required  about  6.5  years  of  CPU  time  on  the  cluster.  The  program  gOpcnMol 
was  used  for  trajectory  visualization.  “  Detailed  analysis  of  individual  trajectories  and  statistical  analysis 
of  the  trajectory  ensemble  was  done  with  programs  specifically  written  for  this  purpose.*  °  available 
from  the  corresponding  author  upon  request. 

One  issue  with  using  quasi-classical  trajectory  methods  is  that  vibrational  energy  is  not  quantized  in 
molecules.  Lack  of  quantization  presumably  has  an  effect  on  how  energy  is  distributed  between 
vibrational  modes  during  trajectories  and  upon  dissociation.10  *  It  is  possible  to  have  trajectories  where 
the  product  vibrational  energy  is  below  the  zero  point  level.  Such  unphysical  trajectories  were  found  to 
be  negligible  in  our  trajectories,  presumably  because  we  are  looking  at  relatively  high  internal  energies 
(sec  fable  2).  reducing  the  cnor  associated  with  treating  the  motion  classically. 

As  mentioned  above,  temperatures  used  for  simulations  are  much  higher  than  the  detonation 
temperature,  which  was  -1500  K  estimated  from  the  temperature  of  the  Bunsen  burner*  used  in 
experiments*0  (For  reference,  the  experimentally  measured  DNB  decomposition  temperature  is  around 
405  K).1  fable  2  lists  the  estimated  vibrational  and  rotational  energies  at  various  temperatures.  As  a 
result,  some  trajectories  may  find  decomposition  paths  with  barriers  too  high  to  be  relevant  to 
decomposition  at  realistic  temperatures.  The  trajectory  observables  of  greatest  interest  for  our  purpose 
arc  decomposition  mechanisms.  The  point  is  that  with  enough  trajectories,  all  important  reaction  paths 
will  be  found,  including  those  that  control  decomposition  under  the  conditions  of  interest.  Once 
decomposition  paths  are  identified,  the  corresponding  reaction  energetics  including  activation  barriers  can 
be  calculated,  and  only  those  which  arc  able  to  contribute  need  be  included  in  decomposition  kinetics 
modeling. 

III.  Trajectory  Results 

Before  we  discuss  the  quasi-classical  trajectory  simulation  results,  it  is  useful  to  review  the  nature  and 
time  scales  of  the  trajectories  we  are  simulating.  Trajectories  were  run  at  the  temperature  range  of  40(H) 
to  6000  K.  At  the  temperature  of 4000  K.  10%  of  trajectories  maintained  the  intact  molecular  structure 


when  trajectories  were  terminated  at  4  pscc.  The  remaining  90%  of  trajectories  at  4000  K  and  all 
trajectories  at  5000  and  6000  K  underwent  thermal  decomposition  via  direct,  concerted  or  sequential 
mechanisms  before  the  termination  of  trajectories.  For  most  of  the  decomposition  trajectories,  the  time  at 
which  the  DNB  molecule  begins  decomposition  is  between  0.2  and  3  pscc  -  longer  than  the  periods  of 
most  DNB  vibrational  modes.  The  decomposition  trajectories  can  be  grouped  into  six  classes  as 
summarized  in  fable  3.  Table  3  also  compares  their  product  energetics,  and  relative  efficiencies  at 
different  temperatures.  The  structures  of  DNB  and  all  decomposition  products  observed  in  trajectories 
arc  depicted  in  Figure  I.  All  these  structures  were  optimized  at  the  B3LYP/6-31-*-K}**  level  of  theory. 
Video  files  showing  representative  trajectories  for  each  of  six  classes  arc  prov  ided  in  the  supporting 
information. 

A.  Decomposition  path  1  (production  of  isocyanic  acid  and  NO:) 

By  far  the  most  common  outcome  of  DN  B  trajectories  is  path  I .  which  accounts  34  -  43%  of  all 
trajectories  at  4000  -  6000K.  Path  I  starts  with  elimination  of  NO:  from  DNB.  followed  by  consecutive 
elimination  of  two  isocyanic  acid  (HNCO)  molecules  from  the  derivative  intermediate.  The 
endothermic ity  of  NO:  elimination  is  1.73  cV.  and  that  of  subsequent  elimination  of  two  HNCO 
molecules  is  2.43  cV.  both  of  which  were  calculated  w  ith  respect  to  the  DNB  molecule.  Figure  2 
demonstrates  a  trajectory  (calculated  at  4000  K)  representative  of  path  1 .  The  figure  show  s  the  change  of 
potential  energy  (PH)  during  the  trajectory.  The  oscillations  in  the  PH  reflect  the  vibration  of  the  DNB 
molecule  and  the  products,  including  ZPH.  The  time  scale  of  decomposition  depends  on  the  simulation 
temperature,  but  three  numbers  arc  relevant.  The  time  between  trajectory  starts  and  the  onset  of  NO: 
elimination  is  around  150  fscc.  The  time  for  elimination  of  the  first  HNCO  molecule  from  the  remaining 
structure  is  around  410  fscc.  and  the  time  for  elimination  of  the  second  HNCO  is  around  800  fscc.  In 
addition  to  the  multi-step  elimination  process  indicated  in  Figure  2.  a  small  fraction  of  path  I  trajectories 
(<  10%)  experienced  simultaneous  elimination  of  NO:  and  HNCO. 

B.  Decomposition  path  2  (production  of  HNNOiH) 


Another  common  class  of  decomposition  tra  jectories  belongs  to  path  2,  accounting  for  7  -  26%  of  all 
trajectories  at  4000  -  6000K.  In  path  2,  HNNO2H  is  produced  via  intra-mo  Iccular  hydrogen  atom  transfer 
followed  by  breaking  an  amide  bond.  HNNO:H  may  further  decompose  to  NNO  and  water.  While  the 
first  step  is  endothermic  by  0.99  cV.  the  overall  decomposition  process,  i.c.,  DNB  0:NNHC0NC0  + 
N:()  +  water,  is  exothermic  by  0.93  eV  with  respect  to  DNB.  A  trajectory  representative  of  path  2  is 
illustrated  in  Figure  3,  indicating  the  energy  barriers  associated  with  two  activated  decomposition  steps  of 
which  the  second  one  is  much  more  significant.  Note  that  in  principle  another  pathway  may  produce  the 
same  products  as  path  2.  i.c.,  DNB  completely  decomposes  to  0:NNHC0NHC0  and  HNNO:, 
subsequently  undergoes  inter-molecular  hydrogen  atom  transfer  between  two  fragments  to  produce 
HNNO>H.  However,  this  decomposition  route  involves  breaking  a  CN  bond  with  an  energy  barrier  in 
excess  of  3.55  cV  calculated  at  B3LYP/6-3H-KJ**.  Therefore,  this  decomposition  path  is  less  likely  to 
contribute  to  a  significant  extent  in  typical  detonation. 

C.  Decomposition  path  3  (production  of  nitramidc) 

Path  3  contributes  to  10%  of  all  trajectories  at  4000  K.  and  the  contribution  decreases  at  high 
temperatures  (6%  at  5000  K  and  4%  at  6000K).  Figure  4  shows  a  typical  trajectory  for  path  3.  In  this 
decomposition  path.  DNB  rotates  a  -NHNO:  unit  out  of  the  molecular  plane,  and  then  undergoes  an  intra¬ 
molecular  SN/-reaction  to  form  nitramidc  H^NNO,  by  transferring  a  hydrogen  atom  from  the  central  N 
atom  to  the  neighboring  amine  N,  followed  by  elimination  of  nitramidc.  This  decomposition  path  has 
been  proposed  in  the  literature, 1  “  and  the  product  nitramidc  was  identified  in  the  mass  spectrum  of  DNB 
pyrolysis  products.1  Similar  to  HNNO:H.  nitramidc  may  also  decompose  to  N:0  and  water  in  the  gas 
phase.1 

D.  Other  decomposition  paths 

The  remaining  trajectories  can  be  characterized  as  paths  4,  5,  or  6,  respectively.  These  paths  arc  much 

less  common,  but  they  arc  similar  in  that  all  involve  sequential  decomposition  reactions  and  the  first  step 

is  NO^  elimination.  In  path  4.  an  mtcr-molecular  hydrogen  atom  transfer  reaction  occurs  between  two 

fragment  compounds  0:NN HCONHCON H  and  NO:,  producing  0:NNHCONHNCO  and  nitrous  acid 
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MONO.  As  visualized  in  path  4  trajectories,  when  the  radical  fragments  O^NNHCONHCONH  and  NO: 
(produced  from  a  harricrlcss  dissociation  path,  vide  infra)  separate  -3  A.  roaming  reorientation  becomes 
feasible  and  this  leads  to  intcrmolccular  M  abstraction  by  NO*  from  ()*NNHCONHCONH.  giving 
molecular  product  NO2H  with  large  vibrational  excitation  of  the  nascent  O-H  bond.  'Hus  decomposition 
path  is  to  some  extent  similar  to  a  roaming  radical  reaction  discovered  recently/1  In  some  path  4 
trajectories,  nitrous  acid  inter-converts  to  UNO*.  At  all  temperatures,  only  10%  (or  less)  of  all  trajectories 
lead  to  path  4.  In  path  5.  0*NNHC0NHC0NH  undergoes  extensive  dissociation  to  NO*,  isocyanic  acid, 
and  a  thrce-mcmbcred  ring  intermediate  cyclic-HNCONH.  The  relative  efficiency  of  path  5  increases 
significantly  at  high  simulation  temperatures,  i.e.,  8%  at  40U0  K  and  increasing  to  13%  at  50U0  K  and 
27%  at  60(H)  K.  presumably  because  this  is  the  most  endothermic  reaction  route.  Finally*  oxygen  atom 
transfer  between  ():NNHCONHCONH  and  NO:  was  observed  in  path  6.  leading  to  net  "NO  elimination” 
from  DNB.  However,  path  6  only  occurs  at  5000  K  with  an  efficiency  of  5%,  and  was  not  observed  at 
4000  and  6000  K. 

These  trajectory  results  lead  us  to  predict  that  paths  I.  2  and  3  may  represent  typical  decomposition 
pathways  of  DNB.  particularly  in  the  view  of  the  initial  decomposition  steps.  These  paths  arc  chosen  for 
several  considerations.  First,  paths  1.  2  and  3  are  energetically  more  favorable  and  overall  attained  much 
higher  probabilities  than  others  in  trajectory  simulations.  Second,  the  products  identified  in  trajectories  of 
paths  1-3.  including  HNCO,  N:().  NO:  and  water,  were  all  observed  experimentally.  4  Therefore*  these 
trajectories  arc  in  reasonable  agreement  with  experiments,  and  can  be  used  as  a  guide  to  extract  additional 
mechanistic  insight  with  confidence.  On  the  other  hand.  NO  (one  of  path  6  products)  accounted  for  only 
small  amount  of  pyrolysis  products,  and  HNO:  (one  of  path  4  products)  was  not  observed  in  pyrolysis.1 
Therefore,  contributions  of  paths  4  and  6  are  expected  to  be  negligible.  Path  5.  regardless  of  any  possible 
energy  bamers.  has  an  endothermic ity  of  3.02  cV  and  is  therefore  of  less  importance  as  well. 

The  following  electronic  structure  calculations  and  statistical  analysis  focus  on  the  initial 
decomposition  steps  of  DNB.  i.e..  elimination  of  NO:,  HNNO:  and  H2NNO2.  It  is  worth  noting  that  paths 
4.  5  and  6  arc  all  related  to  path  1  in  that  these  paths  share  a  common  rate-limiting  step,  i.e.,  NO* 


elimination.  Consequently,  all  possible  decomposition  paths  identified  in  trajectory  simulations  arc 
included  in  the  present  statistical  modeling. 

IV.  KKKM  Modeling 
A.  Reaction  coordinate 

DFT  calculation  results  for  the  reaction  coordinate  of  the  initial  decomposition  paths  arc  summarized 
in  Figure  5  with  the  DNB  molecule  shown  at  zero  energy.  All  energetics  are  calculated  at  the  B3LYP*'6- 
31+-+G**  level  of  theory.  The  details  of  the  geometries  for  the  reactant,  intermediate  complex,  TSs  and 
products  arc  available  by  request  to  the  corresponding  author.  Structures  of  various  TSs  were  optimized 
using  a  number  of  likely  high  potential  energy  geometries  in  reactive  trajectories  as  starting  geometries. 
NO:  elimination  in  path  I  is  dynamically  simple,  in  that  no  complicated  rearrangement  is  needed  for  N-N 
bond  scission.  Efforts  were  made  to  locate  the  transition  state  connecting  DNB  to  products 
O2NNHCONHCONH  +  NO:.  However,  various  starting  geometries  corresponding  to  a  dissociating  N-N 
bond  converged,  instead,  to  the  DNB  structure  at  B3 L Y P/6-3 1++G**  and  MP2/cc-pVDZ  levels  of  theory*. 
We  also  performed  a  relaxed  potential  energy  scan  running  along  the  dissociating  N-N  bond  using  the 
MP2'cc-pVDZ  method.  Ihc  N-N  bond  length  (rss)  was  continuously  varied  from  1 .4  to  3. 1  A,  and  all 
coordinates  other  than  r^*  were  optimized  at  each  point.  The  asymptotic  energy  of  the  resulting  PES  scan 
is  nearly  equal  to  product  energy  ofOjNNHCONHCONH  +  NO:,  suggesting  there  is  no  substantial 
reverse  barrier  associated  with  decomposition  to  0>NN HCON HCON H  +  NO:.  It  is  therefore  expected 
that  NO:  elimination  opens  at  the  thermodynamic  threshold  energy  of  1.73  eV  and  to  be  facile.  Note  that 
previous  ah  initio  calculations  using  G2  and  QCISD(Ty/MP2/6-3l  1G**  methods  indicate  that  H2N-NO2 
may  dissociate  to  NH:  and  NO:  without  a  reverse  barrier/* 4  and  our  calculations  are  consistent  with  that 
finding. 

The  only  low  energy  route  for  path  2  located  computationally  appears  to  be  DNB  —»  TS  2 1  — > 
0:NNHC0NC0  +  HNNO2H.  The  reaction  is  endothermic  by  0.99  eV;  however,  the  decomposition  path 
beais  an  energy  barrier  of  1 . 16  eV  with  respect  to  DNB.  Path  3  is  slightly  more  convoluted:  DNB  — ► 
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TS  31  twisted  DNB  ->  TS  32  ->  0,NNHC0NC0  +  H:NNO:.  i.c.f  twisting  the  nearly  planar  DNB 
molecule  via  an  activation  harrier  TS  31  (0.34  eV  above  DNB),  followed  by  intra-molccular  11-atom 
transfer  and  elimination  of  H:NNO:  via  the  second  activation  harrier  I  S  32  (2.13  eV  above  DNB).  One 
issue  lor  interpretation  of  paths  2  and  3  products  is  that  HNNO>H  might  isomenze  to  M:NNO:,  which 
could  affect  the  subsequent  decomposition  products.  Ihc  I  S  22  for  this  isomerization  is  1.18  eV  above 
tire  HNNO:H  energy.  It  is  certainly  possible  that  inter-conversion  might  occur  during  the  decomposition 
at  detonation  temperature.  No  other  low  energy  pathways  were  found  leading  to  decomposition  paths  l.  2 
and  3,  although  we  cannot  exclude  the  possibility  that  such  pathways  cxisL 
2.  KKkM  results 

To  evaluate  whether  the  initial  decomposition  paths  identified  in  trajectory  simulations  can  account 
for  realistic  situations,  the  RRKM  program  was  used  to  calculate  the  rates  for  each  individual  path  at 
difterent  temperatures.  To  calculate  the  RRKM  unimolecular  rates,  sets  of  vibrational  frequencies  and 
rotational  constants  appropriate  for  the  reactant  and  the  transition  states  leading  to  dissociation  are 
required.  These  arc  derived  from  the  frequencies  and  rotational  constants  calculated  at  B3LYP/6- 
31  —  G**.  with  the  frequencies  scaled  by  a  factor  of  0.955  as  suggested  by  recent  work/  The  rotation 
quantum  number  K  was  treated  as  active  in  evaluating  the  rates  A  (/:.  J)  so  that  all  {2J*\ )  /i- levels  are 
counted.*4  i.c., 

h  Z'^NlE-ErV.K)  | 

where  il  is  the  reaction  path  degeneracy.  O  is  the  transition  state  sum  of  states.  .V  is  the  reactant  density  of 
states.  £*o  is  the  unimolecular  dissociation  threshold,  and  Er  and  Ej  arc  the  rotational  energies  for  the 
reactant  and  the  transition  state,  respectively. 

Since  no  reverse  activation  barrier  is  located  for  decomposition  of  DNB  to  LNNNMC'ONHCONH  +• 
NO>.  vibrational  frequencies  appropriate  to  the  transition  state  for  NO:  elimination,  hereafter  termed 
IS  1,  have  to  be  assumed  for  the  RRKM  calculation  of  path  1  decomposition  rates.  Two  limiting  sets  of 
vibrational  frequencies  have  been  chosen  for  TS  I.  4f"4  ‘  At  one  extreme,  the  vibrational  frequencies  of 
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I  S  1  arc  equal  to  those  of  DNB  with  removal  of  only  one  mode  that  corresponds  to  the  reaction 
coordinate  (i.c..  the  stretching  mode  of  the  dissociating  N-N  bond).  As  long  as  there  is  no  reverse 
activation  barrier  for  decomposition  of  DNB  to  0:NNHC()NHCONH  and  NO:,  this  choice  reflects  the 
"tightest*  TS  that  is  possible,  and  would  provide  lower  limits  to  the  unimolccular  rates  for  path  1 .  At  the 
other  extreme,  a  TS  estimate  was  made  which  consists  of  33  frequencies  of  02 NN HCONHCONH ,  3 
frequencies  of  NO:,  and  5  transitional  frequencies.  The  frequencies  of  the  six  transitional  modes  lost 
upon  dissociation  arc  chosen  from  DNB.  One  N-N  stretching  frequency  is  chosen  as  the  reaction 
coordinate  and  removed.  The  five  remaining  transitional  frequencies  (corresponding  to  rotations,  bends 
and  torsions),  which  become  rotations  and  translations  of  completely  dissociated  products,  arc  sealed  by  a 
factor  of  0.5  to  reflect  the  looseness  of  the  transition  state  which  changes  the  entropy  of  activation  (i.c.  the 
difference  in  entropy  betw  een  the  reactant  and  the  TS).  The  choice  of  our  scaling  factor  is  empirical  and 
was  based  on  other's  w  ork.  A  similar  scaling  factor  was  used  by  Armentrout's  group  to  assign  the 
transitional  modes  in  dissociation  of  energized  metal  ion-ligand  complex  ML*.4M  Dissociations  of  these 
ML*  complexes  have  no  reverse  activation  barriers,  and  the  ion-induced  dipole  and  ion-dipole  forces 
result  in  long  range  metal  ligand  interactions.  Iherefore.  the  corresponding  TSs  must  be  truly  "loose" 
TSs.  The  fact  that  the  similarly  scaled  translational  frequencies  were  successfully  used  to  extract 
dissociation  energies  of  ML*  by  Armentrout  el  al ,4<"4  suggests  that  this  scaling  factor  yields  appropriate 
"loose"  TSs  corresponding  to  weak  associations  of  the  products.  In  the  present  work,  we  employed 
"tight"  and  "loose"  TS  1.  respectively,  in  two  sets  of  RRKM  calculations  for  path  1  NO:  elimination. 

The  tightness  of  the  transition  state  significantly  affects  the  reaction  rate  and  thus  the  branching  ratio,  and 
the  resulting  two  sets  of  RRKM  rates  provide  an  estimate  of  the  range  of  path  1  decomposition  rates. 

fable  IV  summarized  the  RRKM  results,  including  the  densities  of  states  (DOS)  for  the  equilibrium 
and  twisted  DNB  structures,  unimolccular  rates  and  decomposition  branching  ratios.  In  the  fable,  kj,  Ay, 
ksi  represent  the  rates  for  crossing  TS  I,  TS  21.  'f S  3 1  leading  from  the  equilibrium  DNB  structure, 
while  A.j,  and  k#  are  the  rates  for  crossing  TS_3 1  and  TS  32  from  the  twisted  DNB  structure. 
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respectively.  The  table  lists  two  values  for  k)  at  each  temperature  (one  in  parentheses),  which  were 
obtained  using  "tight*  and  ""loose"  TS  I .  respectively.  According  to  the  values  of  ksi  and  k.uy  the  inter- 
convention  between  equilibrium  and  twisted  structures  of  DNB  occurs  on  a  time  scale  of  pscc.  which  is 
significantly  shorter  compared  to  decomposition  time  scale.  Therefore,  the  relative  ratio  of  equilibrium 
and  twisted  DNB  structures  could  be  estimated  using  DOS  in  these  structures.5"  As  indicated  in  Table  4, 
the  values  of  DOS  arc  temperature  dependent.  Because  the  energy  of  the  twisted  structure  increases  by 
only  0.26  cV  compared  to  the  equilibrium  one.  the  DOS  ratios  for  cquilibriunv.twistcd  structure  arc  nearly 
1 : 1  from  750K  to  2000  K.  Consequently,  at  this  temperature  range.  DNB  is  best  considered  to  be  a  mix 
of  two  different  structures,  and  the  average  lifetime  of  DNB  is  around  10**  s  at  1000  K.  decreasing  to  10’9 
s  at  1500  K.  Under  this  circumstance,  the  branching  for  DNB  decomposition  can  he  estimated  by  the 
ratio  of  kj9  k>  and  Aj*. 

RRKM  calculations  predict  that,  at  temperatures  below  1500  K.  path  1  (NO:  elimination)  and  path  3 
(nitramidc  elimination)  arc  negligible  due  to  the  high  and  or  tight  transition  states,  and  only  path  2 
(elimination  of  HNNOdl)  would  be  expected.  Paths  1  and  3  begin  to  emerge  at  1500  K  (i.c..  the  typical 
denotation  temperature40),  and  become  significant  at  high  temperature.  Depending  on  the  nature  of  the 
transition  state  adopted  for  the  calculations,  path  1  may  account  for  1  -  9%  of  initial  decomposition  at 
1750  K  and  2  -  1X%  at  2000  K.  For  the  reason  discussed  above,  this  branching  includes  contributions 
from  paths  4.  5  and  6  in  addition  to  path  1 .  llie  branching  ratio  of  path  3  is  estimated  to  be  1%  at  1500  K. 
increasing  to  5%  at  1750  K  and  12  -  15%  at  2000  K.  Since  TS  32  (for  the  formation  of  nitramidc  from 
twisted  DNB)  and  TS  22  (for  the  formation  of  nitramidc  from  isomerization  of  HNNQ>H)  have  a  similar 
bamcr  height,  it  is  anticipated  that  IKNNO*  may  be  produced  via  both  routes  above  1500K. 

V.  Conclusions 

Quasi-c  lassie  al.  direct  dynamics  trajectories  for  thermal  decomposition  of  DNB  were  analyzed  to 
probe  the  decomposition  mechanisms  and  dynamics.  The  trajectories,  calculated  at  the  B3L  Y  P/6-3 1G* 
level  of  theory,  arc  able  to  reproduce  the  DNB  decomposition  products  observed  in  pyrolysis  experiments 


and  reveal  various  decomposition  paths  leading  to  these  products.  RRKM  theory’  was  used  to  verify  the 
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significance  of  these  decomposition  paths  and  calculate  the  branching  ratios  at  real  detonation 
temperatures,  based  on  the  reaction  coordinate  constructed  using  trajectory  results.  According  to  RRKM 
analysis,  elimination  of  HNNOiH.  H:NNO:  and  NO>,  can  be  important  at  typical  detonation 
temperatures.  While  HNNO:H  elimination  clearly  dominates  at  the  temperature  range  of  750  -  20U0  K. 
the  contributions  of  other  two  paths  increase  substantially  with  increasing  temperature.  At  2000K, 
elimination  of  nitramide  can  account  for  more  than  12%  of  DNB  decomposition.  Contribution  of  NO: 
elimination  depends  on  the  "tightness"  of  the  associated  transition  state,  and  would  account  for  up  to  18% 
of  DNB  decomposition  should  the  reaction  occur  via  a  "loose"  transition  state.  By  combining  direct 
dynamics  simulations  and  statistical  modeling,  we  arc  able  to  accrtain  the  nature  of  the  reaction 
mechanisms  involved  in  the  initial  decomposition  of  DNB.  This  work  demonstrates  the  utility  of  direct 
dynamics  trajectory  simulations  as  a  powerful  guide  for  investigating  decomposition  mechanisms  and 
dynamics  of  complex  energetic  molecules,  locating  transition  states  for  many-atom  systems,  and 
constructing  the  appropriate  kinetics  model;  subsequently,  it  is  necessary  to  use  RRKM  or  other  statistical 
theories  to  get  meaningful  results  as  far  as  for  the  branching  ratios  at  experimental  temperatures  arc 
concerned. 
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Table  1.  Calculated  enthalpy  of  detonation  (A,H)  of  DNB  (kJ'mol) 


Methods 

B3LYP/6-3IG 

B3LYP6-3 IG* 

B3LYP/6-3I+G* 

B3LYP/6-31+4G** 

MP2/cc-pVTZ  ‘ 

Adi 

-606 

-605 

-739 

-757 

-863 

'  Ref  2 
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Table  2.  Calculated  DNB  vibrational  and  rotational  energies  at  various  temperatures 


T/K 

298 

500 

750 

1000 

1250 

1500 

1750 

2000 

4000 

5000 

6000 

1 WeV 

02b 

0.65 

1.27 

1.99 

2.76 

3.57 

4.40 

5.24 

12.25 

15.82 

19.41 

Em/  cV 

0.04 

0.06 

0.10 

0.13 

0.16 

0.19 

0.23 

0.26 

0.52 

0.65 

0.78 
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Table  3.  Direct  dynamics  trajectory  simulation  results  of  DNB  decomposition 


Path 

Decomposition  products 

AH' 

Branching  ratio  %  h 

feV 

4000K 

5000K 

6000K 

r 

s-°  H  v° 

1  1  1 

8  8 

*V° 

-Vy< 

o  o 

1.73 

34 

39 

43 

I  / 

v  •  \  - 

>  a 

2.43 

•  some  trajectories  dissociate  to  NO;.  HN-O.-  and  HNCO 
simultaneously 

2 

*v°  H  s/* 

1  1  ' 

°V°  .V 

.'‘vx ' 

"  1  o 

A  .  1! 

H  H  f| 

N 

0.99 

-0.93 

26 

13 

7 

3 

‘V0  „  °y°  °y° 

-VvV 

8  o  8  8  i 

0.26 

10 

6 

4 

■v0  4  °y° 

/^v\  «/N\ 

8  ^ 

0.54 

4 

v-°  „  %/■ 

.yYY 

-  Y  J  .  <W> 
A^*v\  * 

n  j! 

°  i 

1.73 

4 

8 

10 

0  j 

Ti  --v 

B"  Y  \  i 

S  -  o^o 

A 

l.ll 

22 


Table  3.  (Continued) 


Path 

Decomposition  products 

AH' 

Branching  ratio  % 

leW 

4000K 

5000K  6000K 

5 

°v/° 

i  ■?  i  —  i  v  *  v 

./sA'N, 

S  a  8  8 

1.73 

8 

13 

26 

V° 

* 

S 

C  II 

‘  /  \  *  jj 

NH— *«  I 

3.02 

6 

H  °n/ 

AAA- 

a  s 

I  M  ov 

r  '  N  *  V"  — ^ 

H  Y  Y s" 

1.81 

0 

5 

0 

v  H  „ 

1  |  II  -  MO 

m/S^n>^  s„ 

&  8 

0.98 

Sum  of  oilier  dissociation  paths 

0 

16 

10 

Non-reactive 

0 

10 

0 

■ 

“  Values  of  All  were  calculated  relative  to  DNB.  at  the  B3LYP/6-3H-K***  level  of  theory. 


Branching  ratios  were  calculated  based  on  50  trajectories  at  each  temperature. 
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Tabic  4.  RRKM  results  for  initial  decomposition  of  DNB 


Temp !  K 

750 

1000 

1250 

1500 

1750 

2000 

density  of  states  ( i  VcnV1) 

DNB 

1.95  x  10“ 

2.34  x  10" 

5.98  x  10s5 

4.94  x  10s* 

1.71  x  10" 

2.92  x  10“ 

twisted  DNB 

1.90  x  10“ 

2.33  x  10" 

6.04  x  10“ 

5.04  x  10s* 

1.75  x  10" 

3.02  x  10“ 

unitnolecular  rates  (s’* )  * 

k, 

0 

0.6 

1.8  x  10' 

2.1  x  10" 

3.9  x  10’ 

2.9  x  10* 

(0) 

(3.6) 

(1.8  x  I04) 

(2.7  x  10') 

(5.8  x  10*) 

(4.7  x  lO’’) 

k; 

8.5 

9.7  x  10- 

7.9  x  10’ 

1.0  x  10v 

5.5  x  10' 

1.8  x  lO10 

k„ 

3.4  x  ltf’ 

1.5  x  lO1" 

3.5  x  I01" 

6.0  x  lO10 

8.8  x  10'" 

1.2  x  10“ 

kit 

2.4  x  10" 

3 ’x  10" 

3.8  x  10" 

4.4  x  10“ 

4.8  x  10" 

5.1  x  10“ 

k,- 

0 

1.5  x  10  s 

2.1  x  10* 

7.9  x  10" 

2.9  x  10* 

3.4  x  I0V 

branching  ratios  (%)  * 

path  1 

0 

0 

0 

0.2 

0.7 

1.3 

(0.3) 

(2.6) 

(9.1) 

(18.0) 

path  2 

100 

100 

100 

99.0 

94.3 

83.3 

(99.7) 

(96.7) 

(86.3) 

(69.2) 

path  3 

0 

0 

0 

0.8 

5.0 

15.4 

(0.8) 

(4.5) 

(12.8) 

“  ki,  h,  kn  arc  the  rate  constants  for  crossing  TS_ I .  TS_2 1 .  TS_3 1  from  DNB.  and  k.u  and  ku  arc  the 


rate  constants  for  crossing  TS_3 1  and  TS  32  from  twisted  DNB.  respectively.  The  ki  values  in 
parentheses  were  obtained  using  a  "loose"  TS  I . 
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Figure  Captions 

Figure  1  B3LYP/6-31++G**  optimized  structures  of  DNB  and  products  identified  in  trajectories, 

figure  2  A  representative  trajectory  of  decomposition  path  1,  simulated  at  4000  K  using  the  B3LYP/6- 


31G*  method.  Plot  shows  variation  of  potential  energy  during  the  tra  jectory,  and  the  snap 
shots  show  three  sequential  elimination  steps. 

Figure  3  A  representative  trajectory  of  decomposition  path  2.  simulated  at  4000  K  using  the  B3LYP/6- 
3 1G*  method.  Plot  shows  variation  of  potential  energy  during  the  tra  jectory,  and  the  snap 
shots  show  elimination  of  HNNO2H  and  its  subsequent  dissociation  to  N:0  and  water. 

Figure  4  A  representative  trajectory  of  decomposition  path  3,  simulated  at  4000  K  using  the  B3LYP/6- 
3 1G*  method.  Plot  shows  variation  of  potential  energy  during  the  tra  jectory,  and  the  snap 
shots  show  elimination  ofH^NNOi  via  twisting  the  molecular  plane  followed  by  intra¬ 
molecular  H-atom  transfer. 

Figure  5  Schematic  reaction  coordinate  for  the  initial  decomposition  of  DNB.  Energetics  of 

complexes,  TSs,  and  products  arc  derived  from  B3LYP/6-3I++G**  calculations,  including 
ZPE.  For  TSs,  vibrational  modes  corresponding  to  the  imaginary  frequencies  arc  indicated 
by  displacement  vectors. 
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Fig.  2 
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Fig  .4 
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